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Axion Like Particles (ALPs) are predicted to couple with photons in the presence of magnetic 
fields. This effect may lead to a significant change in the observed spectra of gamma-ray sources 
such as AGNs. Here we carry out a detailed study that for the first time simultaneously considers in 
the same framework both the photon/axion mixing that takes place in the gamma-ray source and 
that one expected to occur in the intergalactic magnetic fields. An efhcient photon/axion mixing in 
the source always means an attenuation in the photon flux, whereas the mixing in the intergalactic 
medium may result in a decrement and/or enhancement of the photon flux, depending on the 
distance of the source and the energy considered. Interestingly, we find that decreasing the value of 
the intergalactic magnetic field strength, which decreases the probability for photon/axion mixing, 
could result in an increase of the expected photon fiux at Earth if the source is far enough. We also 
find a 30% attenuation in the intensity spectrum of distant sources, which occurs at an energy that 
only depends on the properties of the ALPs and the intensity of the intergalactic magnetic field, 
and thus independent of the AGN source being observed. Moreover, we show that this mechanism 
can easily explain recent puzzles in the spectra of distant gamma-ray sources, like the possible 
detection of TeV photons from 3C 66A (a source located at z=0.444) by MAGIC and VERITAS, 
which should not happen according to conventional models of photon propagation over cosmological 
distances. Another puzzle is the recent published lower limit to the EBL intensity at 3.6 fj,m (which 
is almost twice larger as the previous one), which implies very hard spectra for some detected TeV 
gamma-ray sources located at z=0. 1-0.2. The consequences that come from this work are testable 
with the current generation of gamma-ray instruments, namely Fermi (formerly known as GLAST) 
and imaging atmospheric Cherenkov telescopes like CANGAROO, HESS, MAGIC and VERITAS. 

PACS numbers: 95.35.+d; 95.55.Ka; 95.85.Pw; 98.70.Vc; 98.70.Rz; 14.80.Mz 



I. INTRODUCTION 



The existence of axions is predicted by the Peccei- 
Quinn mechanism, which is currently the most com- 
pelling explanation to solve the CP problem in QCD 
Moreover, amongst all the valid candidates pro- 
posed to constitute a portion or the totality of the non- 
barionic cold dark matter content predicted to exist in 
the Universe, hypothetical non-thermal axions, or in a 
more generic way, Axion-Like Particles (ALPs), where 
the mass and the coupling constant are not related to 
each other, may represent a good option: they might ex- 
ist in sufficient quantities to account for the estimated 
dark matter density and they might interact very weakly 
with the rest of the particles 2]. There is an additional 
property of ALPs that makes them even more attractive 
and that could have important implications for its de- 
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tectability, i.e. they can oscillate into photons and vice- 
versa in the presence of an electric or magnetic field [3, ,3] . 
This is analogous to that predicted to occur between neu- 
trinos of different flavors, and a similar behavior is ex- 
pected in the case of the recently proposed chameleons 
as well Q. This characteristic is the main vehicle used 
at present to carry out an exhaustive search of ALPs by 
experiments like CAST fi^, PVLAS 0] and ADMX [1]. 

The oscillation of photons to ALPs (and vice-versa) 
could have important implications for astronomical ob- 
servations. This argument was first investigated in the 
optical band by Ref . Q , where authors proposed the ex- 
istence of axions to be the cause of the observed super- 
nova la dimming. In this context, the observed dimming 
might be explained as a result of an efficient photon to 
axion conversion instead of a cosmic acceleration (albeit 
this proposal was rejected some time later due to some 
chromatic problems, pointed out e.g. in Ref. [13] )■ Pho- 
ton/axion oscillations were also studied by the same au- 
thors in Ref. [llj as an alternative explanation for those 
photons arriving Earth from very distant sources at en- 
ergies above the GZK cutoff. 
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Recently, it has been proposed that, if ALPs exist, they 
could distort the spectra of gamma-ray sources, such as 
Active Galactic Nuclei (AGNs) or galactic 

sources in the TeV range , and that their effect may be 
detected by current gamma-ray experiments. In [23| . for 
example, it is stated that also the scatter in AGN lumi- 
nosity relations could be used to search for ALPs. Other 
astrophysical environments have been proposed in order 
to detect ALPs, such as the magnetic field of the Sun [23 |. 
pulsars , the galactic halo [2^ or GRBs and QSOs by 
carefully studying their polarized gamma-ray emissions 
[13, HI]. In particular, these predictions are very relevant 
for gamma-ray astronomy, where recent instrumentation 
developments in the last few years have increased the ob- 
servational capabilities by more than one order of mag- 
nitude. On the ground, we have the new generation of 
Imaging Atmospheric Cherenkov Telescopes (lACTs) like 
MAGIC m, HESS VERITAS [H or CANGAROO- 
III [l^l, covering energies in the range 0.1-20 TeV. In 
space we have Fermi (previously called GLAST) [sl], in 
operation since Summer 2008 and covering energies in 
the range 0.02-300 GeV[74]. 

In this work we revisit the photon/axion mixing, for 
the first time handling under the same consistent frame- 
work the mixing that takes place inside or near the 
gamma-ray sources together with that one expected to 
occur in intergalactic magnetic field (IGMF). In the liter- 
ature, both effects have been considered separately. De- 
pending on the source dimension, magnetic field, ALP 
mass and coupling constant, both effects might produce 
significant spectral distortions, or one effect could be 
more important than the other. In any case, we believe 
that both effects could be relevant and hence need to 
be considered simultaneously. We neglect, however, the 
mixing that may happen inside the Milky Way due to 
galactic magnetic fields. At present, a concise modeling 
of this effect is still very dependent on the largely un- 
known morphology of the magnetic field in the galaxy. 
Furthermore, in the most idealistic/optimistic case, this 
effect would produce an enhancement of the photon flux 
arriving at Earth of about 3% of the initial photon flux 
emitted by the source [l^ . This is in contrast with what 
we found for the IGMFs: although there is also little in- 
formation on the strength and morphology of the IGMFs, 
the derived photon/axion mixing in this case we show to 
be crucial for a correct interpretation of the observed flux. 
It is worth mentioning that we will come to this conclu- 
sion using a conservative value of B=0.1 nG for the IGMF 
strength, well below the current upper limits of ~1 nG. 
We also carry out a detailed analysis of the mixing when 
varying IGMF strength and source distance. We find re- 
sults that differ from previously published ones, and we 
make predictions of effects that have not been noted in 
the literature so far. 

At energies larger than 10 GeV, and especially above 
100 GeV, it will be necessary to properly account for 
the Extragalactic Background Light (EBL) in our IGMF 
mixing calculations. The EBL introduces an attenua- 



tion in the photon flux due to e~e^ pair production that 
comes from the interaction of the gamma-ray source pho- 
tons with infrared and optical-UV background photons 
[3^. Amongst all the EBL models that exist in the lit- 
erature, in this work we will make use of the Primack 
[Hi and Kneiske best-fit ^ EBL models. They repre- 
sent respectively one of the most transparent and one of 
the most opaque models for gamma-rays, but still within 
the limits imposed by the observations. The EBL model 
will play a crucial role in our formalism and results: as we 
will see, the more attenuating the EBL model considered, 
the more relevant the effect of photon/axion oscillations 
in the IGMF. 

We also explore in this work the detection prospects for 
current gamma-ray instruments (Fermi and lACTs). We 
will show that the signatures of photon/axion oscillations 
may be observationally detectable provided light ALPs 
with masses smaller than a given value for typical values 
of the IGMF. In order to study the detection prospects, 
we will propose an observational strategy. We can an- 
ticipate here that the main challenge for our proposed 
formalism to be testable comes from the lack of knowl- 
edge of the intrinsic source spectrum and EBL density. 
However, we note that there is the possibility that we 
could be already detecting the first hints of axions with 
current experiments. In this context, the potential detec- 
tion of TeV photons from very distant (z 0.4) sources 
[H, 113, lis [in , or some works claiming energy spectral 
indices harder than 1.5 for relatively distant (z=0.1-0.2) 
AGNs (43 , already put in a tight spot the conventional 
interpretation of the observed gamma-ray data. As we 
will show, both effects could be explained by oscillations 
of photon into light ALPs using realistic values for the 
involved parameters. 

The work is organized as follows. In Section |TT] we 
describe in detail the photon/axion mixing in both the 
surroundings of gamma-ray sources and in the intergalac- 
tic medium (IGM). Section Hill is devoted to present the 
results obtained when including both mixings under the 
same framework and after considering realistic param- 
eters for well-known AGNs. In Section IIVI we present 
an observational strategy to search for ALPs using the 
most sensitive gamma-ray instruments, namely Fermi 
and lACTs like MAGIC or HESS. Finally, we give our 
conclusions in Section |Vl 



II. THE FORMALISM 

At present, the Peccci-Quinn mechanism remains as 
the most convincing solution to solve the CP violation of 
QCD. As early as in 1978, Weinberg [ll] and Wilczek 
|44| realized independently that a consequence of this 
mechanism is the existence of a pseudo-scalar boson, the 
axion. One generic property of axions is a two-photon 
interaction of the form: 
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FIG. 1: Sketch of the formahsm used in this work, where both mixing inside the source and mixing in the IGMF are considered 
under the same consistent framework. Photon to axion oscillations (or vice-versa) are represented by a crooked line, while the 
symbols 7 and a mean gamma-ray photons and axions respectively. This diagram collects the main physical scenarios that we 
might identify inside our formalism. Each of them are squematically represented by a line that goes from the source to the 
Earth. 



A. Mixing inside and near the source 
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where a is the axion field, M is the inverse of the pho- 
ton/ axion coupUng strength, F is the electromagnetic 
field-strength tensor, F its dual, E the electric field, 
and B the magnetic field. The axion has the important 
feature that its mass ma and coupling constant are in- 
versely related to each other. There are, however, other 
predicted states where this relation does not hold; such 
states are known as Axion Like Particles (ALPs). An 
important and intriguing consequence of Eq. ([T]) is that 
ALPs oscillate into photons and vice-versa in the pres- 
ence of an electric or magnetic field. In fact this effect 
represents the keystone in ongoing ALP searches carried 
out by current experiments. 

In this work, we will make use of the photon/axion 
mixing as well, but this time by means of astrophysical 
magnetic fields. As already mentioned, we will account 
for the mixing that takes place inside or near the gamma- 
ray sources together with that one expected to occur in 
the IGMFs. We will do it under the same consistent 
framework. Furthermore, it is important to remark that 
it will be necessary to include the EBL in our formal- 
ism, in particular in the equations that describe the in- 
tergalactic mixing. Its main effect we should remember 
is an attenuation of the photon flux, especially at ener- 
gies above 100 GeV. We show in Fig. [1] a diagram that 
outlines our formalism. Very squematically, the diagram 
shows the travel of a photon from the source to the Earth 
in a scenario with axions. In the same Figure, we show 
the main physical cases that one could identify inside our 
formalism: mixing in both the source and the IGMF, 
mixing in only one of these environments, the effect of 
the EBL, axion to photon reconversions in the IGMF, 
etc. A quantitative description of the photon/axion mix- 
ing phenomenon in both the source and the IGMFs can 
be found in the next two subsections. 



It has been recently proposed that an efficient conver- 
sion from photons to ALPs (and vice-versa) could take 
place in or near some astrophysical objects that should 
host a strong magnetic field [12]. 

Given a domain of length s, where there is a roughly 
constant magnetic field and plasma frequency, the prob- 
ability of a photon of energy to be converted into an 
ALP after traveling through it can be written as p^, : 



Po = (As sf 



sm^{Aosc s/2) 



(2) 



{Aosc s/2)2 
Here Aqsc is the oscillation wave number: 

AL^(AcA./ + A,;-A,)2 + 4A|, (3) 
Ab that gives us an idea of how effective is the mixing. 



Ab = 



2 M 



1.7 X 10-21 Mn cm-\ (4) 



where Bt the magnetic field component along the polar- 
ization vector of the photon and Mn the inverse of the 
coupling constant. 

AcM is the vacuum Cotton- Mouton term, i.e. 
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where Bcr — "ig/e ~ 4.41 x 10^'^ G the critical magnetic 
field strength (e is the electron charge). 
Api is the plasma term: 



^i^'i^li^ 3.5 X 10-2" f , , I cm 

^ 2E Vl03cm-3/ \ E^ J 



(6) 
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where Wpi = ^ 47rane/TOe = 0.37 x 10~^/zeV-\/ nf./cm^'^ 
the plasma frequency, me the electron mass and rig the 
electron density. 

Finally, Aa is the ALP mass term: 



2E'y 



2.5 X 10 



(7) 



Note that in Eqs. (|4H7|) we have introduced the dimen- 
sionless quantities B,nG = G, Mn = Af/lO" 

GeV and m^ey ~ m/10^^ cV. 

Since we expect to have not only one coherence do- 
main but several domains with magnetic fields differ- 
ent from zero and subsequently with a potential pho- 
ton/axion mixing in each of them, we can derive a total 
conversion probability as follows: 



1 



[1 - exp(-3ArPo/2)] 



(8) 



where Pq is given by Eq. ^ and N represents the number 
of domains. Note that in the limit where N Pq oo, the 
total probability saturates to 1/3, i.e. one third of the 
photons will convert into ALPs. 

It is useful here to rewrite Eq. ([2]) following Ref. [l2[, 
i.e. 



Pn = 



1 



1 + {E,„t/E^)^ 



sm 
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so that we can define a characteristic energy, Ecru , given 
by: 



m? M 
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or in more convenient units: 
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where the subindices refer again to dimensionless quan- 
tities: m^ey = m/neV, Mn = M/IO" GeV and Bq = 
B/ Gauss; m is the effective ALP mass = — w^jl. 
Recent results from the CAST experiment ^ give a value 
of Mil > 0.114 for axion mass nia < 0.02 eV. Although 
there are other limits derived with other methods or ex- 
periments, the CAST bound is the most general and 
stringent limit in the range 10^^^ eV ^ nia ^ 10^^ 
eV. 

At energies below Ecru the conversion probability is 
small, which means that the mixing will be small. There- 
fore we must focus our detection efforts at energies above 
this Ecrit, where the mixing is expected to be large 
{strong mixing regime). As pointed out in Ref. [l^ . in the 
case of using typical parameters for an AGN in Eq. (fTTj) . 
Ecrit will lie in the GeV range given an ALP mass of the 
order of fj,eV. 



To illustrate how the photon/axion mixing inside the 
source works, we show in Figure [2] an example for an 
AGN modeled by the parameters listed in Table Hi] (our 
fiducial model, see Section HIT)) . The only difference is the 
use of an ALP mass of 1 /icV instead of the value that 
appear in that Table, so that we obtain a critical energy 
that lie in the GeV energy range; we get Ecrit = 0.19 
GeV according to Eq. pT|) . Note that the main effect 
just above this critical energy is an attenuation in the to- 
tal expected intensity of the source. However, note also 
that the attenuation starts to decrease at higher energies 
(>10 GeV) gradually. The reason for this behavior is 
the crucial role of the Cotton-Mouton term at those high 
energies, which makes the efficiency of the source mixing 
to decrease as the energy increases (see Eq. and how 
it affects to Eq. Indeed, the photon attenuation in- 
duced by the mixing in the source completely dissapears 
at energies above around 200 GeV in this particular ex- 
ample. On the other hand, one can see in Figure [5] a 
sinusoidal behavior just below the critical energy as well 
as just below the energy at which the source mixing dis- 
sapears due to the Cotton-Mouton term. However, it 
must be noted that a) the oscillation effects are small; b) 
these oscillations only occur when using photons polar- 
ized in one direction while, in reality, the photon fluxes 
are expected to be rather non-polarized; and c) the above 
given expressions are approximations and actually only 
their asymptotic behavior should be taken as exact and 
well described by the formulae. Therefore, the chances 
of observing sinusoidally-varying energy spectra in as- 
trophysical source, due to photon/axion oscillations, are 
essentially zero. 




FIG. 2: Example of photon/axion oscillations inside the 
source or vicinity, and its effect on the source intensity (solid 
line), which was normalized to 1 in the Figure. We used the 
parameters given in Table [IT] to model the AGN source, but 
we adopted an ALP mass of 1 neV. This gives Ecrit = 0.19 
GeV. The dot-dashed line represents the maximum (theoret- 
ical) attenuation given by Eq. ((8)l, and equal to 1/3. 
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B. Mixing in the IGMFs 

The strength of the Intergalactic Magnetic Fields 
(IGMFs) is expected to be many orders of magnitude 
weaker (~nG) than that of the source and its surround- 
ings ('^G). Consequently, as described by Eq. (fTTjl . the 
energy at which photon/axion oscillation occurs in the 
IGM is many orders of magnitude larger than that at 
which oscillation can occur in the source and its vicinity. 
Despite the low magnetic field B, the photon/axion os- 
cillation can take place due to the large distances, since 
the important quantity defining the probability for this 
conversion is the product B x s, as described by Eq 
Assuming B '--^O.l nG (see below), and Mn = 0.114 (co- 
incident with the upper limit reported by GAST [6|), 
then the effect can be observationally detectable {Ecrit < 
1 TeV) only if the ALP mass is m^ < 6 x 10"^'' eV. If 
the axion mass m^ was larger than this value, then the 
consequences of this oscillation could not be probed with 
the current gene ration of lACTs, that observe up to few 
75| . In our fiducial model (see Table ITU) we 
eV, which implies Ecrit = 28.5 GeV. 



tens of TeV 
used ma = 10 
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It is important to stress that at energies larger than 10 
GeV, and especially larger than 100 GeV, besides the os- 
cillation to ALPs, the photons should also be affected by 
the diffuse radiation from the Extragalactic Background 
Light (EBL). The EBL introduces an attenuation in the 
photon flux due to e~e~^ pair production that comes from 
the interaction of the gamma-ray source photons with in- 
frared and optical-UV background photons for the ener- 
gies under consideration [3a |. Therefore, it will be neces- 
sary to modify the above equations to properly account 
for the EBL in our calculations. These equations can be 
found in Ref. [ll[, where the photon/axion mixing in the 
IGMF was also studied, although for other purposes and 
a different energy range. We note that the same equa- 
tions were also used in Ref. |l3| to study for the first 
time the photon/axion mixing in the presence of IGMFs 
for the same energy range that we are considering in this 
work. 

There is little information on the strength and mor- 
phology of the IGMFs. As for the morphology, several 



authors reported that space should be divided into sev- 
eral domains, each of them with a size for which the mag- 
netic field is coherent. Different domains will have ran- 
domly changing directions of B field of about the same 
strength J^, J6J . The IGMF strength is constrained to 
be smaller than 1 nG [47|, which is somewhat supported 
by the estimates of ~0.3-0.9 nG that can be inferred [i^ 
from recent observations of the Pierre Auger Observa- 
tory [4^ . On the other hand, there is controversy on the 
possibility of generating such a strong magnetic field. De- 
tailed simulations yield IGMFs of the order of 0.01 nG 
so that they can later reproduce the measured B fields in 
nearby galaxy clusters |50l . Isij . Given this controversy, 
we decided to use a mid- value of 0.1 nG in our fiducial 
model (Table HD. 

In our model, we assume that the photon beam prop- 
agates over N domains of a given length, the modulus of 
the magnetic field B roughly constant in each of them. 
We will take, however, randomly chosen orientations, 
which in practice will be also equivalent to a variation in 
the strength of the component of the magnetic field in- 
volved in the photon/axion mixing. If the photon beam 
is propagating along the y axis, the oscillation will occur 
with magnetic fields in the x and z directions since the 
polarization of the photon can only be along those axis. 
Therefore, we can describe the beam state by the vector 
{lx,^z, a)- The transfer equation will be [Tl| : 



7, I = e^'^v e^^y + Ti e^^'' + T2 e^'^ ] 7. I (12) 
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e being the angle between the x-axis and B in each single domain. S a dimcnsionless parameter equal to: 

, BX^ / B \/10"GeV\/A^ 

= — -— — 0.11 ' ' ' ' ' 
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FIG. 3: Effect of intergalactic photon/axion mixing on photon and ALP intensities versus distance to the source, computed for 
our fiducial model, i.e. for 3C 279 and those parameters given in Table HIl but taking B= 1 uG, and using the Primack EBL 
model. The black thick solid line represents the total photon intensity, while the blue dotted line is the ALP intensity. The 
photon intensity as given only by the EBL (i.e. without including photon/axion mixing) is shown as the red dashed line. Top 
panels: mixing computed for Mn = 4 GeV and an initial photon energy of 50 GeV (left), 500 GeV (middle) and 2 TeV (right); 
bottom panels: Mn = 0.7 GeV and same energies than top panels. 



that represents the number of photon/axion oscillations 
within the mean free path of the photon Notice that 
if there was no EBL, the quanta beam would be equipar- 
titioned between the ALP component and the two photon 
polarizations after crossing a large number of domains. 
However, the EBL introduces an energy dependent mean 
free path A-y for the photon. 

Amongst all the EBL models that exist in the liter- 
ature, we chose the Primack [s^ and Kneiske best-fit 
[sj to fix the parameter. They represent respectively 
one of the most transparent and one of the most opaque 
models for gamma-rays, but still within the limits im- 
posed by the observations (galaxy counts for the lower 
limit and observations of distant blazars for the upper 
one). The model proposed by Kneiske et al. was initially 
disfavored by some TeV observations of distant AGNs, 
using the assumption that the intrinsic spectral index 
needs to be softer than 1.5 (see (E^H^). On the other 
hand, in the literature we also find work where this as- 
sumption is strongly criticized, as reported by [s^ . Issj 
and especially in [42]. Therefore, we will consider the 
Kneiske best-fit EBL model as still valid. In the case of 
the Primack EBL model, we did not take into account 
the redshift evolution of the EBL, which effect may be 
particularly important for sources located at z > 0.3. Al- 



though one of the objects we take as the reference, 3C 
279, is located at a larger redshift, the difference in the 
final photon intensity when using the Primack EBL with 
or without redshift evolution is still small (below 15% for 
the relevant energies). This is not true for the Kneiske 
EBL model, for which the differences might be as high as 
70% for some of the energies under consideration. There- 
fore, it became necessary to account for such effect in this 
case. 

From each EBL model, we obtain A^ as the distance 
given by the so-called gamma-ray horizon for the energy 
considered. Additionally, we have to take into account 
that the energy of each photon will change continuously 
for a photon traveling towards us from cosmological dis- 
tances, due to the cosmological redshift. This effect may 
have a very important role in the calculations of the pho- 
ton/axion mixing, since e.g. for a source at a distance of 
1000 Mpc (i.e. z ~ 0.3) every photon arrives at Earth 
with 30% less energy. We account here for this effect 
for the first time by computing at each step (distance) 
the new energy of the photon due to cosmological red- 
shift, and then using this new energy as the input energy 
needed for the calculation of A.^. We did not include 
in the formalism, however, those secondary photons that 
may arise from the interaction of the primary source pho- 
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tons with the EBL. 

To illustrate how the mixing in the IGM works, we 
show in Figure [3] various examples of the evolution of the 
total photon and ALP intensities as a function of the dis- 
tance to the source when varying some of the critical pa- 
rameters, using the Primack EBL model in all cases. We 
use the parameters listed in Table |TT1 that corresponds 
to our fiducial model, but using an IGMF strength of 
1 nG (instead of 0.1 nG), which is still consistent with 
upper limits. For a photon with an initial energy of 50 
GeV (left top panel) and coupling constant Afn = 4, 
which yields Ficrit ~ 100 GeV, there is not significant 
photon/axion oscillations. Since the role of the EBL is 
almost negligible at these energies, the total photon in- 
tensity remains almost constant traveling to the Earth. 
For 500 GeV photons (middle top panel), the total pho- 
ton intensity initially decreases as expected due to the 
EBL absorption, but also an early extra attenuation due 
to a photon to ALP conversion is clearly observed. At the 
same time, the ALP intensity, which was initially equal 
to zero (we neglect here the mixing inside the source for 
simplicity), grows rapidly. At larger distances the ten- 
dency in the total photon intensity is just the opposite; 
the intensity increases slightly, since an efficient ALP to 
photon reconversion (although operative since the very 
beginning) is taking place and becomes relevant specially 
at these large distances, where the expected photon in- 
tensity is already very low due to the EBL absorption. 
In the case of photons with higher initial energy (e.g. 2 
TeV, right top panel in Fig. [3]), the expected attenua- 
tion due to the EBL becomes very important even for 
small distances from the source, which makes more rele- 
vant the impact of ALP to photon reconversions on the 
photon intensity. As a result, the photon/axion mixing 
implies an enhancement in the photon intensity at almost 
all distances. The situation changes when using a slightly 
higher coupling constant, but still within the CAST con- 
straints (sec bottom panels in Fig. [3]). In this case, both 
the attenuation and the enhancement in intensity become 
more pronounced, as expected. For relatively small dis- 
tances, the photon/axion mixing produces an attenuation 
in the photon flux, while for relatively large distances, the 
mixing produces an enhancement in the photon flux. The 
same argument is essentially valid for any initial photon 
energy, the results only changing depending on the rela- 
tive relevance of the EBL in each case, which will modify 
the distance at which the photon-intensity enhancement 
starts occurring. 



III. RESULTS 

Up to now, previous works have focused only in study- 
ing the photon/axion mixing either inside the source or 
in the IGMFs. Instead, for the flrst time we carried out 
a detailed study of the mixing in both regimes under the 
same consistent framework. We neglect for the moment, 
however, the mixing that may happen inside the Milky 



Way due to its galactic magnetic flelds. We believe that 
a concise modeling of this effect is still very dependent 
on the largely unknown morphology of the B field in our 
Galaxy. In the most idealistic/optimistic case, in which 
30% of the photons convert to ALPs within the source 
and 10% of the ALPs convert to photons in the Milky 
Way, as reported in Ref. 15], this effect would produce 
an enhancement of the photon flux arriving at Earth of 
about 3% of the initial photon flux emitted by the source. 

As mentioned in the previous section, in order for the 
photon/axion oscillation to be observationally noticeable 
by current instruments, that is Ecrit < 1 TeV, we need 
ALP masses smaller than 6 • 10~^'^eV for typical val- 
ues of the IGM. Larger ALP mass values translate into 
higher Ecrit (e.g. m^ = 10~^ eV would yield Ecrit ^ PeV 
in the IGM, when using B '--^O.l nG) making the effect 
of this oscillation undetectable by the current gamma- 
ray instruments (Fermi and lACTs). In scenarios with 
heavy ALPs, then the only effect detectable would be 
an attenuation caused by photon/axion oscillation in the 
source, which would be of about 30% in the most op- 
timistic case [l2|. It is worth mentioning that the po- 
tential photon/axion oscillation in the Milky Way could 
produce measurable effects only for ALP masses smaller 
than 10~* eV, as mentioned in Ref. p^ . 

We use an ALP mass of 10"^" eV in our fiducial model 
(Table nil), which implies Ecrit ^ 30 GeV in the IGM (for 
B^O.l nG) and Ecrit ~ 1 eV within the source and its 
vicinity (B~l G). Consequently, both effects need to be 
taken into account. 



A. Photon/axion oscillation in our framework 

In this section we show the results obtained when tak- 
ing into account the mixing inside the source and in the 
IGMF simultaneously. Since we expect the intergalactic 
mixing to be more important for larger distances, due to 
the more prominent role of the EBL, we chose two dis- 
tant astrophysical sources (as our benchmark AGNs) that 
are relatively well characterized at gamma-ray energies; 
namely the radio quasar 3C 279 (z=0.536), most distant 
detected gamma-ray source at the VHE range, and the 
BL Lac PKS 2155-304 at z=0.117. In order to compute 
the photon/axion oscillation in the source we used the 
parameters reported in Ref. [5^ for 3C 279 and Ref. [s^l 
for PKS 2155-304. As for the size of the region with B 
field (where photons can convert to ALPs) we chose a 
region 10 times larger than the radius of the gamma-ray 
emitting blob given in the above mentioned references, 
the reason being that the blob radius represents only a 
lower limit for the region where B is confined. We note 
that this parameter, as well as the number of domains 
where the B is coherent, play an important role in the 
photon attenuations due to the photon/axion mixing in 
the source. Table U shows the different attenuations that 
are obtained when varying the size of the region where B 
is confined (what we called "B region") and the lengths 
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of coherent B domains inside that region. One can see 
that, once the number of domains is fixed, the photon 
attenuation increases when increasing the size of the "B 
region" . On the other hand, when fixing the size of the 
"B region" and scanning the size of the domains we find 
that, as we increase the number of domains, the attenua- 
tion increases until the size of the domain is "too small" . 
At this point, the probability of photon/axion conversion 
is almost zero for the single domains, which reduces the 
overall photon/axion conversion. 

TABLE I: Maximum attenuations due to photon/axion oscil- 
lations in the source obtained for different sizes of the region 
where the magnetic field is confined ("B region") and differ- 
ent lengths for the coherent domains. Only length domains 
smaller than the size of the B region are possible. The B 
field strength used is 1.5 G (see Table HHl . The photon flux 
intensity without ALPs was normalized to 1. In bold face, is 
the attenuation given by our flducial model. 



B region (pc) 


Length domains (pc) 


3x10-* 3x10-" 0.03 0.3 


0.3 
0.03 

3x10-^ 


0.84 0.67 0.67 0.75 
0.98 0.84 0.77 - 
0.99 0.98 



We summarize in Table |TT] the parameters we have con- 
sidered in order to calculate the total photon/axion con- 
version in both the source (for the two benchmark AGNs) 
and in the IGM. As we already mentioned, these values 
represent our fiducial model. 

TABLE IL Parameters used to calculate the total pho- 
ton/axion conversion in both the source (for the two AGNs 
considered, 3c279 and PKS 2155-304) and in the IGM. The 
values related to 3C 279 were obtained from Ref. [5^, while 
those ones for PKS 2155-304 were obtained from Ref. [STj . 
As for the IGM, ed,int was obtained from [s^, and Bi„t was 
chosen to be well below the upper limit typically given in the 
literature (see discussion in the text). This Table represents 
our fiducial model. 





Parameter 


3C 279 


PKS 2155-304 




B (G) 


1.5 


0.1 


Source 


ed (cm ^) 


25 


160 


parameters 


L domains (pc) 


0.003 


3 X 10"" 




B region (pc) 


0.03 


0.003 




z 


0.536 


0.117 


Intergalactic 


ed.int (cm~^) 


10-^ 


10-^ 


parameters 


B,„t (nG) 


0.1 


0.1 




L domains (Mpc) 


1 


1 


ALP 


M (GeV) 


1.14 X 10'" 


1.14 X 10"' 


parameters 


ALP mass (eV) 


10-"' 


10-"' 



The effect of existence of ALPs on the total photon 
flux coming from 3C 279 and from PKS 2155-304 (using 
the fiducial model presented in Table HI]) can be seen in 
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FIG. 4: Effect of photon/axion conversions both inside the 
source and in the IGM on the total photon flux coming from 
3C 279 (z=0.536) and PKS 2155-304 (z=0.117) for two EBL 
models: Kneiske best-fit (dashed line) and Primack (solid 
line). The expected photon fiux without including ALPs is 
also shown for comparison (dotted line for Kneiske best-flt 
and dot-dashed line for Primack). 

Figure SI We carried out the calculations for the two 
EBL models cited above: Kneiske best-fit and Primack. 
The inferred critical energies for the mixing in the source 
are Ecrit = 4.6 eV for 3C 279 and Ecrit = 69 eV for 
PKS 2155-304, while for the mixing in the IGM we ob- 
tain Ecrit = 28.5 GeV. The photon attenuation due to 
photon/axion mixing inside the source is 16% for 3C 279 
and 1% for PKS 2155-304, as can be seen above their re- 
spective critical energies in Figure 21 On the other hand, 
the photon attenuation due to photon/axion oscillation 
in the IGM is 30% for the distance of both sources, and 
it occurs at the same critical energy. The role of the EBL 
is negligible at this low energy (i.e. below ^100 GeV), 
which means that the intensity curves for the two EBL 
models agree to this energy. 

The situation changes above '^100 GeV, where the 
photon attenuation due to the EBL is noticeable. At 
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this point, the results depend substantially on the source 
distance and the EBL model used. A stronger photon 
attenuation is obtained for the Kneiske best-fit model 
against the Primack EBL model, as expected. Because 
the strong photon attenuation due to the EBL, the ALPs 
that later convert to photons imply a further enhance- 
ment of the expected photon flux. Therefore, as one can 
notice from Figure IH the existence of ALPs translates 
into a relatively small (~30%) intensity attenuation at 
low energies and a large intensity enhancement (several 
orders of magnitude, depending on the energy range, dis- 
tance of the source and chosen EBL model) at high en- 
ergies. 

In order to quantitatively study the effect of ALPs on 
the total photon intensity, we plot in Figure [5] the dif- 
ference between the predicted arriving photon intensity 
without including ALPs and that one obtained when in- 
cluding the photon/axion oscillations (called here the ax- 
ion boost factor). Again, this was done for our fiducial 
model (Table and for the two EBL models described 
above. The plots show differences in the axion boost fac- 
tors obtained for 3C 279 and PKS 2155-304 due mostly 
to the redshift difference. 

In the case of 3C 279, the axion boost is an atten- 
uation of about 16% below the critical energy (due to 
mixing inside the source). Above this critical energy and 
below 200-300 GeV, where the EBL attenuation is still 
small, there is an extra attenuation of about 30% (mix- 
ing in the IGMF). Above 200-300 GeV the axion boost 
reaches very high values: at 1 TeV, a factor of ~7 for the 
Primack EBL model and ~340 for the Kneiske best-fit 
model. As already discussed, the more attenuating the 
EBL model considered, the more relevant the effect of 
photon/axion oscillations in the IGMF, since any ALP 
to photon reconversion might substantially enhance the 
intensity arriving at Earth. We note that the axion boost 
factor may vary when changing the parameters we used 
to model the source (as shown in Table [l| and the IGM 
(see next section). The results we find in this work are 
in disagreement with those reported by De Angelis et 
al. [isl. We always find that the photon intensity below 
200-300 GeV decreases when including the oscillation to 
ALPs regardless of the ALP and/or IGM parameters, 
while De Angelis et al. find that the photon intensity 
increases for a large ran ge o f the phase space they tried 
(see their Fig. 1 in Ref. [iSl). At those low energies the 
photon attenuation due to pair conversion in the EBL is 
relatively low (see Fig. [J]) and thus the few ALPs that 
convert to gamma photons do not imply any substantial 
relative increase in the photon intensity. On the other 
hand, 1/3 of the photons oscillate to ALPs, which causes 
a substantial decrement in the amount of gamma pho- 
tons with respect to those we would have in the absence 
of ALPs. Therefore, we think it is very difficult to get a 
photon enhancement at energies ~100 GeV. On the other 
hand, the axion boost factors we find at high energies 
(>300 GeV) are substantially lower than those obtained 
in [l^. As an example, in the case of a Kneiske best-fit 



EBL model with B=l nG, we find a boost ~4 at 500 
GeV, whereas De Angelis et al. obtain ~20 for the same 
photon energy and the same redshift (note that, in order 
to carry out a one-to-one comparison with that work, we 
also used Mu = 4, as they do). One of the reasons for the 
discrepancy in the axion boost factors is the used EBL 
model. We noted that the EBL model shown in Fig. 1 of 
Ref. (l3| is substantially more attenuating than the one 
from Kneiske best-fit EBL model, which is the one we are 
using. Consequently, the axion boost factors reported in 
Ref. [31 are larger than the ones they would have ob- 
tained if they had used the Kneiske best-fit EBL model. 
Besides that, it is not clear to us whether the change in 
photon energy due to cosmological redshift (see Section 
IIIBI) was taken into account in [IjI ; this is not explicitely 
mentioned in their work. 

In the case of PKS 2155-304, the situation is differ- 
ent from that of 3C 279 due to the very low photon- 
attenuation at the source and, mostly, due to the smaller 
source distance. The low redshift location decreases the 
impact of the EBL absorption and thus the effect of the 
relative photon- flux enhancement due to photon/axion 
oscillation. When using B=0.1 nG, the axion boost fac- 
tor is larger than 1 only for the Kneiske best-fit model 
and only above 1.3 TeV. In the case of Primack EBL, the 
axion boost factor is always smaller than 1, thus implying 
no photon- flux enhancement. Note however that the 30% 
drop in the photon intensity occurs at the same energy as 
that of 3C 279. This drop in the photon intensity should 
occur at the same energy for all sources located at rel- 
atively medium redshifts (0.1<z<0.3). For very nearby 
sources (z <0.05), the energy drop should still be the 
same since it only depends on the ALPs properties and 
the strength of the IGMF. However, the magnitude of 
the drop will decrease. This is thus a very distinctive 
and easily testable prediction of this mechanism. We 
will discuss this issue again in Section [TVl 

We want to stress another interesting feature that the 
photon/axion oscillation in the IGMF produces in the 
source spectra at VHE (>100 GeV) of distant sources 
(z>0.1). As one can notice from Fig. O the axion boost 
factor starts increasing at few hundred GeV (when the 
EBL becomes important), and consequently it will make 
the source spectra to look harder than they are actually. 
This happens for both AGNs, yet at slightly different en- 
ergies: 100 GeV for 3C 279 and 300 GeV for PKS 2155- 
304. As shown in Figure. [5l the hardening of the VHE 
spectra occurs for both (very different) EBL models that 
we used, and hence a very robust prediction of this mech- 
anism being at work. Such a hardening of the spectra was 
already predicted in Ref. [lH| for several AGNs located 
at redshifts 0.1-0.2. While in our work the effect is due 
mostly to the photon/axion oscillation in the IGMF, in 
(ist the effect is due to photon/axion oscillation within 
the source (up to 30% attenuation of the photon flux) 
and the one that occurs in the galactic magnetic fields 
of the Milky Way (up to 10% conversion probability). It 
is worth mentioning here, however, that when using the 
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FIG. 5: Boost in intensity due to ALPs for the Kneiske best- 
fit (daslied line) and Primack (solid line) EBL models, com- 
puted using the fiducial model presented in Table HIl for 3C 279 
(z=0.536) and PKS 2155-304 (z=0.117). 



parameters (essentially B strength and size of the "B re- 
gion" ) for the modeling of the gamma-ray emission from 
AGN sources, we find that the attenuation in the source 
due to photon/axion conversion is relatively low; 16% for 
our model of 3C 279 and 1% for that of PKS 2155-304. 
These low photon-flux attenuation (equivalent to ALP- 
enhancement) would decrease significantly the effect of 
the mechanism proposed in 15]. 

Finally, it is worth mentioning that we checked that 
our results are robust against the randomness of the B 
field. We ran 100 different realizations of the same phys- 
ical scenario, randomly varying the orientation of B in 
each coherent domain and each realization. We did so for 
the four cases studied along this work, i.e. 3C 279 and 
PKS 2155-304, Primack and Kneiske best-fit EBL mod- 
els. Furthermore, we repeated the same exercise using 
0.1, 1 and 10 Mpc as the length of the coherent domains 
in order to explore the dependence of our results on this 
parameter. In all cases, we chose our fiducial value of 0.1 



nG for the B field strength. We found that the maximum 
differences are typically well below 10%, implying that 
the results obtained are not sensitive to the randomness 
of the B field. We increased the number of realizations to 
1000 for some cases and found no differences with respect 
to the results obtained with 100 realizations. 

A larger effect on the computed axion-boost factors oc- 
curs when changing the size of the domains being used. 
The computed axion-boost factors are sensitive to choice 
of the size of the coherent domains to be used. Together 
with the choice of the EBL model (which is also uncer- 
tain), the choice of the domain sizes modifies the results 
obtained by factors of a few. 



B. The impact of changing B 

A very interesting result has been found when vary- 
ing the modulus of the intergalactic magnetic field. In 
Ref. [Tsj . the intergalactic photon/axion mixing was re- 
jected arguing that its effect on the final intensity at 
Earth would be negligible when using a more realistic 
value for B, which should be substantially lower than 
the value of 1 nG adopted in Ref. fl3\. However, as was 
shown in the previous section, when using B=0.1 nG we 
find significant effects even for sources located at red- 
shifts as low as Z'^0.1. In order to quantify the impact of 
changing the IGMF strength, we plot in Fig.[B]the result 
of varying B in our fiducial model by one order of mag- 
nitude (above and below). We do that for both 3C 279 
and PKS 2155-304. 

In the case of 3C 279, we see in the left top panel of 
Fig. [5] that higher intensities (or equivalently, higher ax- 
ion boost factors in the right top panel), are obtained 
when using B=0.1 nG instead of taking B=l nG. This 
seems to contradict the intuitive idea of getting higher 
intensities for larger magnetic fields, that make the pho- 
ton/axion mixing more efficient. The reason for this re- 
sult is the strong attenuation due to the EBL. If the 
photon/axion mixing is efficient, then many ALPs con- 
vert to photons which soon disappear due to the EBL 
absorption. Consequently, if the source distance is large, 
we end up having a very small number of photons arriv- 
ing to the Earth. On the other hand, if the photon/axion 
mixing is not that strong, then we can keep a higher num- 
ber of ALPs traveling towards the Earth, which act as a 
potential reservoir of photons. When decreasing B to 
0.01 nG, then the axion boost factors are lower than for 
the other two cases. On the other hand, in the case of 
PKS 2155-304, we see that the highest axion boost fac- 
tors are obtained with B=l nG, because the source is 
not as distant as 3C 279. If we had considered a source 
located at a much further distance than 3C 279, then 
we would have found the highest axion boost factors for 
B=0.01 nG. 

In summary, higher B values do not necessarily trans- 
late into higher photon flux enhancements. There is al- 
ways a B value that maximizes the axion boost factors; 
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FIG. 6: Same as in Figures |4] and [5] but for different values of IGMF. Upper panels: 3C 279 using those parameters listed in 
Table HIl only changing B. Lower panels: Same exercise for PKS 2155-304, using the corresponding parameters that can be 
found in the same Table Ull 



this value is sensitive to the source distance, the consid- 
ered energy and the EBL adopted model. 



C. The impact of using the smallest photon/ ALP 
coupling constant 

The most stringent limits on the ALP-photon coupling 
constant were derived using the non-detection of gamma- 
rays (by the Solar Maximum Mission Gamma-Ray Spec- 
trometer) from the supernova (SN) 1987A during the ~10 
seconds time window defined by the neutrino burst. This 
outstanding event allowed several authors in 1996 to set 
lower limits to the inverse of the coupling constant Mn 
to values larger than 1 and 3 Those limits are 
only valid for ultralight ALPs. In both works the value 
nia < 10~^ eV is quoted, although this value holds only 
for some specific situations. Indeed, a more robust value 
is nia < 10~^^ eV (see Refs. [l3,[ll|): i-e. the energy be- 
low which the exact value of the ALP mass is irrelevant 



because the "plasma frequency" dominates (see defini- 
tion of ALP effective mass in Section III Ap . Various au- 
thors (see Refs. [HI, HI]) used Mn — 4 when dealing with 
TUa < 10^° eV. Since the ALP mass in our fiducial model 
is 10"^" eV, and hence close to this limit, we decided to 
repeat the calculations using this value for A/, which is 
35 times larger than the value we used in the previous 
sections (see Table lll| . 

Before we continue, it is worth pointing out that the 
limits to the ALP-photon coupling constant given in 
Refs. [H, [13] are subject to large uncertainties that are 
not fully discussed in those papers. Both the flux of ALPs 
produced in the SN explosion and the back-conversion 
of ALPs to gamma photons can vary by large factors, 
and hence the upper limits computed with those num- 
bers have to be taken with caveats. 

The calculated flux of ALPs produced and released 
during the SN explosion depends on the knowledge of 
the size, temperature and density of the proto-neutron 
star. Those numbers are subject to large uncertainties 
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because we still do not know how stars explode. Even 
though there is general agreement that the ultimate en- 
ergy source is gravity, the relative roles of neutrinos, 
fluid instabilities, rotation and magnetic fields continue 
to be debated. In particular, back in the 90s it was be- 
lieved that neutrinos would be able to reheat the outgo- 
ing shock-wave and produce the explosion. Nowadays, 
with far more powerful computer simulations, we know 
that neutrino-driven explosions are only possible when 
the star has a small iron core and low density in the sur- 
rounding shells, as being found in stars near or below 10 
solar masses [la]. The progenitor of SN1987A was a blue 
supergiant and hence it is expected to be somewhere be- 
tween 10-50 solar masses. A possibility to explain those 
explosions might require the proper inclusion of rotation 
and magnetic fields (see Refs. [13, HO, IMi and references 
therein). Both B field and rotation are present in stars 
as well as in pulsars, which are the products of successful 
SN explosions; thus it is very natural to consider them 
in SN explosion models. In particular, the rotation of 
the proto-neutron star can change substantially the tem- 
perature and, specially, the density of the inner core; in 
[2l| it is shown that the density can vary by more than 
one order of magnitude, which would change by a similar 
factor the flux of ALPs being produced. Refs. [IB, [I3| 
did not consider such level of complexity (and uncertain- 
ties) in the parameters used to compute the flux of ALPs, 
mostly because 15 years ago we lacked that knowledge. 

On the other hand, the back-conversion of ALPs to 
photons relies on the structure of the galactic magnetic 
field which is, again, not well known. Different models 
predict B fields that could differ substantially and hence 
they would predict different values for the amount of 
gamma photons we would obtain for a given flux of ALPs. 
This is clearly shown in Fig. 1 from Ref. [15], where the 
probability of ALP-photon conversion is given for var- 
ious locations of the sky. Therefore, even if we could 
accurately predict the number of ALPs from SN1987A, 
the number of photons would be subject to large uncer- 
tainties. 

Therefore, we conclude that the limit in the inverse of 
the ALP-photon coupling constant given in Refs. [H, [13] 
is subject to large (orders of magnitude) uncertainties, 
and thus the limit given by the CAST collaboration re- 
mains as the most robust one up to date. However, for 
the sake of comparison with other works, we computed 
the axion-boost factors when using Afn ~ 4 eV. This is 
shown in Figure[7]for both 3C 279 and PKS 2155-304 for 
two values of the B field, 0.1 nG and 1 nG. For this low 
coupling constant, the effect due to the photon/ ALP os- 
cillation in the source is negligible. The effect due to pho- 
ton/ALP oscillation in the IGMF in not negligible, but 
substantially lower than the one shown in the previous 
section. Besides, such effect shows up at larger energies 
now (see Eq. HI]); 100 GeV and 1000 GeV respectively 
for 1 nG and 0.1 nG. 
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FIG. 7: Boost in intensity due to ALPs for the Kneiske best-fit 
and Primack EBL models, computed using the fiducial model 
presented in Table|lI]for 3C 279 (z=0.536) and PKS 2155-304 
(z=0.117), but with Mil = 4 GeV and B=0.1 nG (dashed 
and solid lines for Kneiske best-fit and Primack EBL models 
respectively) and B=l nG (dot-dashed and dotted lines). 



IV. DETECTION PROSPECTS FOR FERMI 
AND lACTS 

As mentioned in Section |TT1 for photon/axion coupling 
constants close to the current published limits, and for 
realistic ALP mass values, the energy at which the pho- 
ton/axion oscillation starts to become important is ex- 
pected to lie in the gamma-ray range. Consequently, 
the combination of the Fermi/LAT instrument and the 
lACTs, which cover 6 decades in energy (from 20 MeV 
to 20 TeV) is very well suited to study the photon/axion 
mixing effect. Because of the rapid change in the pre- 
dicted photon intensity attenuation close to Ecru, the 
energy resolution of the instrument is very relevant in 
order to detect the ALP signatures. The Fermi/LAT in- 
strument has an energy resolution of about 10%, whereas 
lACTs, above 150 GeV, have an energy resolution of 
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about 20-25%. The photon intensity attenuation we 
found in this work goes from to 30% due to the mixing 
in the source, plus essentially 30% due to the mixing in 
the IGM. This implies that one needs to be able to de- 
termine photon fluxes with a precision better than 10%. 
Such level of precision might not be achievable for ener- 
gies >10 GeV with Fermi, or for energies >1 TeV with 
lACTs due to a low photon counting (which depends ob- 
viously on the brightness and hardness of the gamma-ray 
sources). On the other hand, if the source is emitting at 
~TeV energies and is located at large distances, the pho- 
ton/axion oscillation in the IGM could translate into an 
intensity enhancement of more than one order of magni- 
tude (see Figures m El 0, which should be certainly easy 
to detect with current lACTs. 

Therefore, if we accurately knew the intrinsic spectrum 
of the sources and/or the density of the EBL, we should 
be able to observationally detect ALP signatures for a 
wide range of the parameter space (photon/axion cou- 
pling constant and ALP mass). The main problem is 
that we do neither know accurately the intrinsic spec- 
trum of the sources nor the EBL density. Thus the po- 
tential detection of those ALP signatures become quite 
challenging, but not impossible. In order to study this 
scenario, we propose the following strategy: 

1. Observe several AGNs located at different redshifts, 
as well as the same AGN undergoing diS'erent flar- 
ing states (low/high fluxes), at different energy 
ranges, from radio to TeV. This is important be- 
cause the modeling of the gamma-ray emission de- 
pends critically on the emission at lower energies 
(specially infrared, optical, UV and X-rays), and 
also because we do not know a priori the energy 
at which the photon/axion oscillation will start to 
operate. 

2. Try to describe the observational data with 
"conventional" theoretical models for the broad 
band emission (Synchrotron Self-Compton, Exter- 
nal Compton, Proton synchrotron, etc) and the at- 
tenuation of the gamma-rays in the EBL (Primack, 
Kneiske best-fit or other EBL models) . The current 
models (sometimes very simplistic) will definitely 
require some modifications to fit the observational 
data. 

3. Look for intensity drops in the residuals ("best- 
model" -data). We want to stress that the drop in 
the photon fiux due to the attenuation in the IGM 
only depends on the IGMF and the properties of 
the ALPs (mass and coupling constant), i.e. it is 
independent of the gamma-ray sources. Therefore, 
a detection of such photon flux drops at the same 
energy in a numerous of different sources would be 
a clear signature for the existence of photon/axion 
oscillation, because we do not expect that the in- 
trinsic spectrum from different sources (or same 
source at different flux levels) have a rapid drop 



of ~30% in the emission at the same energy. The 
detection (or no-detection) of this photon intensity 
drop imphes a constraint for the product to^ • Mn. 
In this specific search, the Fermi-LAT instrument 
is expected to play a key role since it will detect 
thousands of AGN sources located at various red- 
shifts (up to z^5), and at energies where the photon 
absorption due to the EBL is not important. 

4. Look for intensity enhancements in the residuals 
(i.e. "best-model" -data). This should occur at the 
highest energies (E>300 GeV) and thus only de- 
tectable with lAGTs. The origin of the potential 
photon flux "excess" might be due to a wrong EBL 
model and/or wrong model for the source emis- 
sion, the last being very important because it intro- 
duces differences between the different sources (or 
the same source under different activity levels). In 
this case what we need is to detect distant (z > 0.2) 
sources at the highest possible energies (>1 TeV). 
The current EBL models are already very close to 
the minimum possible photon density limits from 
galaxy counts; that is, we cannot make them much 
more transparent. That implies that the detection 
of TeV photons from a source that is at redshift 0.5 
(like 3G 279) could not be explained with conven- 
tional physics, regardless of the intrinsic spectrum 
of the source. This would be a strong hint for the 
existence of photon/axion oscillation. If the same 
effect was observed for different sources at differ- 
ent redshifts, we could try to parameterize the ef- 
fect by varying the ALP parameters (and/or the 
IGMF strength). If that parameterization could be 
done successfully, then we would not only have a 
very strong hint for the existence of ALPs, but also 
would be able to constrain the available parameter 
space (coupling constant and ALP mass). 

The detection of ALPs is not trivial (and cannot be 
done with just few sources), but it is certainly possible, 
as we have shown above. Along these lines, it is worth 
mentioning that we might be already starting to see hints 
of the existence of ALPs from the gamma-ray spectra 
of cosmological sources. Very recent works already pose 
substantial challenges to the conventional interpretation 
of the observed source spectra from several distant AGN 
sources. On the one hand, the VERITAS GoUaboration 
recently claimed a detection of gamma-rays above 0.1 
TeV (the highest detected energies are not yet reported) 
coming from 3C 66A [? ], an intermediate- frequency- 
peaked BL Lac object located at redshift 0.444. This 
claim coincides with the detection, at GeV energies, of 
this source in active state by the Fermi-LAT instru- 
ment during the same time window 59]. In addition, 
the MAGIC Collaboration reported gamma-rays above 
1 TeV coming from a location consistent with the posi- 
tion of 3G 66A [1^. Those observations would confirm 
the earlier claims from the Crimean Astrophysical Ob- 
servatorys GT-48 lACT of several detections from this 
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source above 0.9 TeV [1^ |4l|. Those detections were 
not confirmed by the HEGRA and Whipple telescopes, 
which are more sensitive instruments, but which observed 
the source at different time windows [60l [Glj . As men- 
tioned above, a detection of TeV photons from a source 
located at z=0.444 would pose serious problems to con- 
ventional models of photon propagation over cosmolog- 
ical distances, where the high energy gammas are ex- 
pected to disappear due to pair electron-positron produc- 
tion in the EBL. On the other hand, the recent published 
lower limits to the EBL at 3.6 microns [i^l, which is al- 
most twice larger as the previous ones, enhances even fur- 
ther the attenuation of gamma-rays at TeV energies and 
thus increases even more the magnitude of the mystery. 
Furthermore, as reported in [i^l, this fact extends the 
problems to sources located at medium redshifts (z=0.1- 

0. 2) whose intrinsic energy spectra appear to be harder 
than previously anticipated. Those observations present 
blazar emission models with the challenge of producing 
extremely hard intrinsic spectra (differential spectral in- 
dex in the spectrum smaller than 1.5) in the sub- TeV to 
multi-TeV regime. As mentioned in the previous section, 
the photon/axion oscillation in the IGM would naturally 
explain these two puzzles; the detection of TeV photons 
from very distant (z~0.5) AGNs, and the apparent hard- 
ening of the spectra for relatively distant (z>0.1) AGNs. 

However, it is worth mentioning that the above re- 
ported puzzles might still be explained with conventional 
physics, as well as uncertainties in the published num- 
bers. The measured redshift of 3C 66A could be wrong 
[H, ill, or the TeV photons reported by MAGIC could 
come from a neighboring source (a radio galax y), 3C 66B, 
which has never been detected in gamma-rays [39| . As for 
the blazars with intrinsic spectra harder than 1.5, there 
is currently quite some controversy. Some authors claim 
that spectra harder than 1.5 could be possible (see e.g. 
[63,[6l,[6l,[73]), while others state that spectra should be 
always softer than 1.5 (see for instance |63l. Isol [tH. [t^). 
and use this value to set upper limits to the EBL density 
at infrared frequencies [ssl . l63l . [t^ . . An argument in 
favor of the latter is the fact that EGRET never mea- 
sured spectra harder than 1.5 at energies below 10 GeV, 
where the EBL does not distort the gamma-ray spectra, 
for any of the almost 100 detected AGNs [zi]. 

Finally, we would like to note that the capabilities of 
detecting the mentioned signatures will increase signifi- 
cantly with the new generation of ground instruments, 

1. e. MAGIC II or HESS II (with lower energy thresholds, 
expected to operate in 2009), CTA and AGIS (with even 
lower energy thresholds and higher sensitivity at multi- 
TeV energies) and HAWC (higher sensitivity at multi- 
TeV energies, large duty cycles). 



V. CONCLUSIONS 

If ALPs exist, then we should expect photon to ALP 
conversions (and vice-versa) in the presence of magnetic 



fields. This photon/axion mixing will occur in gamma- 
ray sources as well as in the IGM. We have explored in 
detail both mixing scenarios together in the same frame- 
work. The main conclusions on this work can be sum- 
marized as follows: 

• If photons oscillate into ALPs in the IGM, then 
photon/axion mixing in the source is also at work 
for lower photon energies. In this picture, both ef- 
fects should be taken into account using the same 
framework, since they will be governed by the same 
set of physical parameters (ALP mass and coupling 
constant). In the case of ALP masses ma >> 10^^" 
eV, the energies at which the photon/axion os- 
cillation occur in the IGMF are >> 1 TeV and 
thus not detectable with current gamma-ray instru- 
ments. In those cases the photon/axion oscillation 
in the source would be the only effect that could 
potentially be detected. 

• The photon/axion oscillation in the source (and 
its vicinity) can produce photon-flux attenuations 
up to 30%, as previously stated in the literature 
[H, llBl- However, when using available models for 
gamma-ray emitting blob regions to set values of 
the B field strength and the size of the region where 
the conversion can take place (we took a radius 10 
times the size of the blob), we obtain photon-flux 
attenuations that are significantly lower. 

• The photon/axion oscillation in the IGM produces 
a photon-flux attenuation up to 30% below the en- 
ergies at which the EBL is important (but above 
Ecrit for the oscillation to be efficient). If the source 
redshift is larger than ^0.1, this drop in intensity 
should be about 30% and it shows up in all sources 
at the same energy. Hence, it presents relatively 
easy signature of the presence of ALPs. The Fermi- 
LAT instrument is expected to play a very impor- 
tant role in this search, since it is expected to detect 
thousands of AGN sources located at various red- 
shifts (up to z=5), and at energies where the EBL 
is not relevant. The detection of such a photon in- 
tensity drop would set the value for the product 

• M, under the assumption of a given IGMF 
strength. If such an intensity drop is not seen in 
the spectra, lower limits could be set. 

• Above energies at which the absorption of gamma- 
rays in the EBL become important, the pho- 
ton/axion oscillation in the IGMF could produce 
both attenuation and enhancement in the photon 
flux, depending on the source distance and energy 
under consideration. 

• We find that decreasing the intensity of the IGMF 
strength does not necessarily decreases the photon- 
flux enhancements (axion boost factors). For a 
source located at z=0.5, B=0.1 nC produces higher 
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photon-flux enhancements that B=l nG. This re- 
sult is somewhat unexpected since stronger B fields 
allow for a more eflicient photon/axion mixing. 
The reason for this result is the strong attenuation 
due to the EBL. If the photon/axion mixing is ef- 
ficient, then many ALPs convert to photons which 
soon disappear due to the EBL absorption. Conse- 
quently, if the source distance is large, it ends up 
having a very small number of photons arriving at 
the Earth. On the other hand, if the photon/axion 
mixing is not that efficient (lower B field), then 
there is a higher number of ALPs traveling (to- 
wards the Earth), which act as a potential reservoir 
of photons. The net balance between the two pro- 
cesses is sensitive to the source distance, the energy 
considered and the EBL intensity. Given those pa- 
rameters, there is always a B value that maximizes 
the photon flux enhancements. 

We have shown that the signatures of photon/axion os- 
cillations may be observationally detectable with current 
gamma-ray instruments (Fermi/LAT and lACTs). Since 
photon/axion mixings in both the source and the IGM 
are expected to be at work over several decades in en- 
ergy, it is clear that a meticulous search for ALPs in the 
(sub)GeV-(multi)TeV regime will be greatly enhanced by 
means of a joint effort of Fermi and current lACTs. 

The main challenge in such detection comes from the 
lack of knowledge in conventional physics; namely the 
intrinsic source spectrum and EBL density and the in- 
tensity and configuration of the intergalactic magnetic 
field. In other words, the effect of the photon/axion os- 
cillations could be attributed to conventional physics in 
the particular source and/or propagation of the gamma- 
rays towards the Earth. However, we believe that such 



photon/axion oscillations could be studied using several 
distant AGNs located at different redshifts, as well as 
the same distant AGN detected at distinct activity levels. 
The signatures of such effect being attenuations (at rela- 
tively low energies) and/or enhancements (at the highest 
energies) in the photon fluxes, that could be visible in 
the residuals from the "Best-Model-Fit" and the obser- 
vational data. 

Recent work, like the potential detection of TeV pho- 
tons from very distant (z ^ 0.4) sources, or those ones re- 
porting energy spectral indices being harder than 1.5 for 
relatively distant (z=0.1-0.2) AGNs, already pose sub- 
stantial challenges to the conventional interpretation of 
the observed gamma-ray data. Both effects could be ex- 
plained by oscillations of photons (using ~0.1 nG for the 
IGMF strength) into light ALPs {nia < 10" i° eV) with 
a photon/axion coupling constant close to current upper 
limits (Mil ^ 0.114). 
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